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Abstract 

Solving numerically the equations of motion for the effective lagrangian 
describing supersymmetric QCD with the SU(2) gauge group, we find a 
menagerie of complex domain wall solutions connecting different chirally asym- 
metric vacua. Some of these solutions are BPS-saturated walls; they ex- 
ist when the mass of the matter fields does not exceed some critical value 
m < m* = 4.67059 .... There are also sphaleron branches (saddle points of 
the energy functional). In the range m* < m < to** ~ 4.83, one of these 
branches becomes a local minimum (which is not a BPS - saturated one). At 
m > m.**, the complex walls disappear altogether and only the walls connect- 
ing a chirally asymmetric vacuum with the chirally symmetric one survive. 
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1 Introduction 



Supersymmetric QCD is the theory involving (in a simplest version) a gauge vector 
supermultiplet V and a couple of chiral matter supermultiplets belonging to the 
fundamental representation of the gauge group SU(N). The lagrangian of the model 
reads 

C = ^-Ltt I d 2 6 W 2 + H.c.j + -/ d 2 9d 2 9 £W 

- (^jd 2 6SfS f + K.c) , (1.1) 

where Sf = £f g S 9 (for further conveniences, we have changed a sign of mass here 
compared to the standard convention). The dynamics of this model is in many 
respects similar to the dynamics of the standard (non-supersymmetric) QCD and, 
on the other hand, supersymmetry allows here to obtain a lot of exact results [Q. 

Like in the standard QCD, the axial Ua(1) symmetry corresponding to the chiral 
rotation of the gluino field and present in the tree-level lagrangian is broken 
by anomaly down to Z 2 n- This discrete chiral symmetry can be further broken 
spontaneously down to Z 2 so that the chiral condensate < Tr {A a A Q } > (X a is a 
two - component Weyl spinor describing the gluino field) is formed. There are N 
different vacua with different phases of the condensate 

< Tr A 2 > = Z e 27rik/N , k = 0, . . . , N - 1 (1.2) 

It was noted recently that on top of N chirally asymmetric vacua (|1.2| ), also a 
chirally symmetric vacuum with zero value of the condensate exists. 

The presence of different degenerate physical vacua in the theory implies the 
existence of domain walls — static field configurations depending only on one spatial 
coordinate (z) which interpolate between one of the vacua at z = — oo and another 
one at z = oo and minimizing the energy functional. As was shown in 0, in many 
cases the energy density of these walls can be found exactly due to the fact that 
the walls present the BPS-saturated states. The key ingredient here is the central 
extension of the N = 1 superalgebra @, [3J] 

{QaQp} = 4 /^V{ 

A domain wall is a configuration where the integral of the full derivative in the RHS 
of Eq.( |1.3| ) is non-zero so that the standard N — 1 SUSY algebra in the wall sector 
is modified. A BPS-saturated wall is a configuration preserving 1/2 of the original 
supersymmetry i.e. a configuration annihilated by the action of two certain real 
linear combinations of the original complex supercharges Q a , Q a 0. 

3 Such BPS-saturated walls were known earlier in 2-dimensional supersymmetric theories (they 
are just solitons there) |6| and were considered also in 4-dimensional theories in stringy context 
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Combining (|1.3|) with the standard SUSY commutator {QayQ 13 } = 2(0^)^^ 
and bearing in mind that the vacuum expectation value of the expression in the 
square brackets in Eq.(|1.3|) is zero due to Konishi anomaly |5) [], it is not difficult to 
show that the energy density of a BPS-saturated wall in SQCD satisfies a relation 
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where the subscript ±00 marks the values of the gluino condensate at spatial infini- 



ties. The RHS of Eq.( |1.4|) presents an absolute lower bound for the energy of any 
field configuration interpolating between different vacua. 

The relation ( |1.4j ) is valid assuming that the wall is BPS-saturated. However, 
whether such a BPS-saturated domain wall exists or not is a non-trivial dynamic 
question which can be answered only in a specific study of a particular theory in 
interest . 

In Refs. j|, [TIJ this question was studied for the SU{2) gauge group in the frame- 
work of the effective low energy lagrangian due to Taylor, Veneziano, and Yankielow- 
icz [IllJ. It is written in terms of the composite colorless chiral superfields 



32tt 



r Tr {W 2 }, X 2 = S f S 



f 



;i.s) 



The effective potential is rigidly fixed from the requirement that the conformal and 
the chiral anomaly of the theory ( |1 . 1|) under consideration are reproduced correctly. 
There is some freedom of choice for the kinetic term which we will discuss later in 
the paper. The effective TVY lagrangian presents a Wess-Zumino model with the 
superpotential 
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The corresponding potential for the lowest components 0, \ °f the superfields $, X 
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( from now on, with few exceptions, we will measure everything in the units of 
A-sqcd = A) has three non-trivial minima: 
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i See recent |9[ for detailed pedagogical explanations. 
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There are also the minima with the inverse sign of x (and the appropriately chosen 
phase of 0), but they are physically the same as the minima ( |1.9|) : the vacuum 
expectation values of the gauge invariant operators < Tr A 2 > = (327r 2 /3) < 3 > 
and < s'sf > = < \ 2 > (s? is the squark field) are the same. The vacua ( |i~9| ) 
involve a non-zero gluino condensate ( |1.2| ). The vacuum ( |1.8| ) is chirally symmetric. 

To study the domain wall configurations, we should add to the potential ( |1.7|) 
the kinetic term which we choose in the simplest possible form C^ n = \d<f)\ 2 + \ dx\ 2 
and solve the equations of motion with appropriate boundary conditions. In Ref . [|J , 
it was done for the walls connecting a chirally asymmetric vacuum with the chirally 
symmetric one. It turned out that the solution exists for any value of m and that it 
is always BPS - saturated, i.e. the walls could be found solving not the equations 
of motion, but the more simple first order BPS equations (we again refer a reader 
to Ref.|| for details) 

d z <f) = ±<9W/<90, d zX = ±dW/dx (1.10) 

The surface energy density of these walls coincides with its lower bound ( |1.4| ). Ex- 
plicitly, in physical units 



e = 4 V / (mA| QCD )/3 (1.11) 

However, there is also a different kind of walls which connect different chirally 
asymmetric vacua in (|1.9|) []. The BPS equations with the corresponding boundary 
conditions were solved numerically in [H|. It was shown that the solutions exist 
only in the limited range of masses m < m* = 4.6705 .... At larger masses, complex 
BPS domain walls do not exist. A kind of phase transition occurs. 

In this paper, we perform a more detailed study of the complex walls. We show 
that on top of the upper BPS branch found in Ref.|3j], also a lower BPS branch 
exists. Here the terms "upper" and "lower" refer not to energy ( the energy of both 
branches is, of course, the same and is given by the BPS bound 2e ), but to some 
other parameter: the value of \<f>\ at the middle of the wall. When the mass is small, 
the trajectories can be found analytically. For the lower branch, the fields in the 
middle of the wall are very small in this case so that the wall presents a very loose 
bound state of two real walls. When the mass rises, |0| m iddie also rises. At m = m», 
two branches fuse together. 

We also studied numerically the wall solutions of the second order equations of 
motion. Some of the solutions are the BPS walls described above, but some are not. 
At m < m*, there are two new sphaleron branches. They present saddle points of the 
energy functional. The upper (in the same sense as above) sphaleron branch presents 

5 We will call these walls complex walls in contrast to the real walls studied in Ref. [Q. For the 
walls separating the symmetric and an asymmetric vacua, the phase of the fields does not change. 
It can be set to zero which considerably simplifies the equations. For the walls separating different 



vacua in (1.9), the trajectories 4>(z) and x( z ) are essentially complex, and the problem is technically 



more intricate. 
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a "mountain pass" between two different BPS valleys. The lower sphaleron branch 
presents a mountain pass between the lower BPS valley and the field configuration 
corresponding to two real BPS walls separated at infinite distance. At m = when 
two BPS valleys fuse together, the upper sphaleron solution coincides with the BPS 
one: a mountain pass disappears and there is only one local minimum of the energy 
functional. This local minimum (and the corresponding solution for the equations 
of motion) exists also at somewhat larger masses, in the range m* < m < m**. It is 
not, however, a BPS - saturated minimum anymore, its energy somewhat exceeding 
the BPS bound. At m = « 4.83, the upper branch fuses together with the lower 
sphaleron branch. At larger masses, no complex domain wall solutions exist. [] 

The plan of the paper is the following. Before proceeding with numerics, we 
discuss in the next section the limits m — > and m — > oo where the problem can be 
studied analytically. In Sect. 3, we describe in details the numerical solutions of the 
BPS equations in the intermediate mass range. In Sect. 4, we solve the second order 
equations of motion and display the existence of two sphaleron branches. The last 
section is devoted to the discussion of physical relevance of our results and to some 
speculative remarks on the issues which are not yet clear. 



2 Exploring the limits. 

The situation is greatly simplified when the mass parameter m in the lagrangian 
( |1 . ip is large or small (compared to A). Consider first the case of large masses. 
In this case, one can integrate the heavy matter fields out and write the effective 
lagrangian for the composite chiral superfield $. Technically, one should use the 
Born-Oppenheimer procedure and to freeze down the matter field x so that the large 
second term in the potential ( |1.7|) disappear. Proceeding in supersymmetric way, 
we get X 2 = 4$ 3 /3m. Substituting it in the first term, we obtain the Veneziano- 
Yankielowicz effective lagrangian [12| which is the Wess-Zumino model for the single 



chiral superfield $ with the superpotential 



4 , 
W = -$ 3 
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(2.1) 



where A| yM = ^3mA 5 S Q CD /A. The expression (|2.1|) and the corresponding expres- 
sion for the potential 

U(<p) = \dW/d<f)\ 2 = 16|0 4 ||ln0 3 | 2 (2.2) 

(here is measured in the units of A$ym) is not yet well defined: the logarithm has 
many sheets, and one should specify first what particular sheet should be taken. An 

6 To avoid misunderstanding, we emphasize that two chirally asymmetric vacua can always be 
connected in two steps, via the symmetric one. That is there is always a "composite" solution 
presenting a combination of two real walls separated at infinite distance. 
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accurate analysis taking into account the fact that the topological charge 

1 

16tt 2 

in the original theory (|1.1|) is quantized to be integer reveals that the true potential 
is glued out of two pieces 0, M : 



v = — I Tr{G 2 G 2 }d 4 x (2.3) 



ln0 3 = ln|0 3 | +iarg(0 3 ), arg(0 3 ) G (-7r/2,7r/2) 
ln0 3 = ln|0 3 | +i[arg(0 3 ) -vr], arg(0 3 ) G (vr/2,3vr/2) (2.4) 

(remind that only the field <p z ~ Tr A 2 , not <fi itself has a direct physical meaning). 
Quite an analogous situation holds in the Schwinger model: the true bosonized 
lagrangian (in the Schwinger model it is just equivalent to the original theory ) is 
glued out of several branches when taking into account the effects due to quantization 
of topological charge . The potential ( |2.2| ) has three degenerate minima. Two of 



them with < Tr A 2 > oc < 3 > = ±1 are chirally asymmetric and the third one 
with < <p > = is chirally symmetric. 

The domain walls in this model were studied in Q. It turned out that only 
the real domain walls interpolating between a chirally asymmetric and the chirally 
symmetric vacua exist in this case. Nontrivial complex domain wall solutions which 
do not go through zero are absent. Speaking of the real walls, they are BPS saturated 
and can be obtained solving the first order equations d<fi = ±dW/d<p. The solution 
is "almost" analytical, it is expressed via integral logarithms. 

Consider now the case of small masses. In this case, chirally asymmetric vacua 
are characterized by large values of the matter scalar field x ■ Again, the theory 
involves two different energy scales, and one can tentatively integrate out heavy 
fields and to write the Wilsonean effective lagrangian describing only light degrees 
of freedom. Proceeding in the Born-Oppenheimer spirit, we should freeze now the 
heavy field in the potential ( |1.7| ) so that the large first term in the potential acquire 
its minimum (the zero) value. In contrast to the large mass situation, this can now 
be achieved in two ways: either by setting <ft = or by setting <f) 3 x 2 = 1. In the first 
case, we will obtain the effective lagrangian describing the dynamics of the chirally 
symmetric phase (it is just the lagrangian of free light chiral field X involving \ and 
its superpartner). 

The second choice results in the lagrangian describing the dynamics of the chi- 
rally asymmetric phases. It is the lagrangian of the Wess-Zumino model with a 
single chiral superfield X and a non-trivial superpotential 

W = 2 —- -X 2 . (2.5) 

3X 2 2 V ; 

This lagrangian is well known and was obtained earlier from instanton and/or from 
holomorphy considerations [1J. The corresponding potential U = \dW/dx\ 2 has 
two different non-trivial minima at 



< x 2 > = ±xt = ±V4/3m (2.6) 
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When m < 1, large expectation value ( |2.6| ) results in breaking down the gauge 
symmetry of the original theory by the Higgs mechanism: the theory is in the Higgs 
phase. 

A domain wall interpolating between two vacua ( |2.6|) is BPS-saturated. The 
solution can be found analytically M 

1+i 4m(z-z ) 

X (z) = x» = 2.7 

y/\ _|_ e 8m(z-z ) 

where z$ is the position of the wall center. We see that \x{ z )\ stays constant and 
only the phase is changed. 

In this approach, we are not able, however, to study the real domain walls 
interpolating between the chirally symmetric and a chirally asymmetric vacua. Such 
a wall corresponds to going through a high energy barrier separating two kind of 
vacua. It is a remarkable consequence of supersymmetry and of the related BPS 
condition (|1.4j) that the energy of such a wall is still not large. 

It was shown in M that in the limit m — > 0, the real wall trajectory can also 
be found analytically. In a funny way, the transition from the chirally symmetric 
to a chirally asymmetric vacuum goes in two stages. On the first stage, only the 
field \x\ is changed rising from zero at z = — oo to its maximal value x* at z = 0, 
x(z) = X* e7nz ' i while the field stays frozen at zero. This first part of the wall 
has the width ~ 1/m and carries the half of the wall total energy. At z > 0, the 
trajectory abruptly turns. x( z ) is n °t changed anymore, but <f>(z) changes until it 
levels off at its value in the asymmetric vacuum state 0* = x* 2 ^ ■ The equations 
for this second stretch have exactly the same form as the BPS equations in the pure 
super symmetric Yang-Mills theory. The second part of the wall carries another half 
of the total energy, but is much thinner than the first part, its width being of order 
m -1 / 6 . It is on this second stage when the potential barrier between symmetric and 
asymmetric vacua is penetrated. 

Thereby, the physical situation in the limits m — > and m — > oo is somewhat 
different. In both cases, we have 3 different vacua. However, an analog of nontrivial 
complex walls (|2.7|) connecting different asymmetric vacua which are present at small 
masses, are absent when the mass is large. In the latter case, only real walls are 
present. It is therefore very interesting to understand what happens in between, at 
intermediate values of masses, and how the transition from one regime to another 
occurs. That was the main motivation for our study. 

Generally, one should study the theory with the potential (11.71 ). The status of this 
effective theory is somewhat more uncertain than that of ( |2.5| ) — for general value 
of mass, the TVY effective lagrangian is not Wilsonean; light and heavy degrees of 
freedom are not nicely separated. But it possesses all the relevant symmetries of the 
original theory and satisfies the anomalous Ward identities for correlators at zero 
momenta. We think that the use of the TVY lagrangian is justified as far as the 
vacuum structure of the theory is concerned. 
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3 Solving BPS equations. 



Let us choose the positive sign in Eq. ( |1.10|) and try to solve it with the boundary 
conditions 



3m N 1 / 6 



-oo 



/ 6m\ 



oo 



/ 4 \'/ 4 
X( " 00) = v3m" 



^ /3 (xf' 



(3-1) 



(the negative sign in ( |1.10 ) would describe the wall going in the opposite direction 
in z). The solution of the equations (|1.10|) with the boundary conditions ( |3.1|) has 
the fixed energy density which due to (|1.4j) is twice as large as the energy density of 
the real wall Thereby a complex BPS wall may be thought of as a marginary 

stable bound state of two real walls. 

Technically, it is convenient to introduce the polar variables \ = pe ia , <ft = Re lf3 - 
Then the system ( |1.10|) (with the positive sign chosen) can be written in the form 



d z p 
d z a 
d z R 



—mp cos(2ct) + ^ cos(3/3) 

msin(2a) — |^-sin(3/3) 

2R 2 [cos(3/3) \n(R 3 p 2 ) - sin(3/3)(3/3 + 2a)} 

-2R [sin(3/3) ln(R 3 p 2 ) + cos(3/3)(3/3 + 2a)} 



(3.2) 



The wall solution should be symmetric with respect to its center. Let us seek for 
the solution centered at z = so that 

p( z ) = p(-z), R(z) = R(-z), a{z) = tt/2 - a{-z), p(z) = -tt/3 - P(-z) (3.3) 



Indeed, one can be easily convinced that the Ansatz ( |3.3| ) goes through the equations 

(ID- 

The system ( |1.1(J ) has one integral of motion |J: 



Indeed, we have 



0, 



Im W(4>, x) — const 



■ -r^o z x 



(3.4) 



dx 



dW 


2 


dW 


d<j> 


+ 


dx 



d z W 



It is convenient to solve the equations ( p.2|) numerically on the half-interval from 
z = to z = oo. The symmetry ( |3.3| ) dictates a(0) = 7r/4, (3(0) = —n/6. Then the 
condition (|3.4j ) [in our case Im W(0, x) — due to the boundary conditions (|3.1|) 1 
implies 



4i?(0) s 



{ln[i?(0) 3 p(0) 2 ] - l} +mp(0) 2 = 



(3.5) 



7 A logical possibility could be that symmetric equations have a couple of asymmetric solutions. 
We searched for such but did not find any, however. 
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m — m * 



Figure 1: Mismatch parameter A as a function of mass 

Thus, only one parameter at z = [say, -R(O)] is left free. We should fit it so that 
the solution would approach the complex minimum in Eq. ( |1.9| ) at z — ► oo. 

It turns out that the solution of this problem exists, but only in the limited 
range of m. If m > m* = 4.67059 . . ., the solution misses the minimum no matter 
what the value of -R(O) is chosen. This is illustrated in Fig. [I] where the "mismatch 
parameter" 

A = nunmin^/|x(;z) - i(4/3my/ 4 \ 2 + \<f>(z) - e~ i,r /3(3 m /4)i/6| 2 ( 3 . 6 ) 

is plotted (in the double logarithmic scale) as a function of mass. The dependence 
A(m) fits nicely the law 

AbpsM = 0.56(m - m,) - 44 (3.7) 

It smells like a critical behavior but, as A is not a physical quantity, we would not 
elaborate this point further. 

At m = 4.67059 or at smaller values of mass, the solution exists, however. 
Moreover, it turns out that, at m < m*, there are two different solutions: the upper 
branch with the larger value of -R(O) and the lower branch with the lower value. In 
Fig. H we plotted the dependence of rj — R(0)/R(oo) on m for both branches. We 
see that, at m = m*, two branches are joined together. This is the reason why no 
solution exists at larger masses. 

Let us discuss now what happens with these two branches in the small mass 
limit. Consider first the upper branch. The profiles of the functions p(z) and R(z) 
for three values of mass: m = 0.01, m = 2.0, and m = 4.67059 for the upper branch 
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in the whole interval — oo < z < oo are plotted in Figs. [3],f|. We see that, for small 
masses, R(z) and p(z) are virtually constant. This is not surprising, of course. For 
small masses, the solution approaches, as it should, the analytic solution (|2.7|) (with 

= x" 2/3 )- 

Linearizing at small nonzero masses the equations ( |3.2| ) with respect to small 
deviations of the solutions Ap(z), Aa(z), etc. from the analytic result (|2.7| ) ( and 
<p = x~ 2 ^ 3 ) valid in the limit m — > 0, it is possible to determine the dependence 
r](m) in the small mass region. The derivation in the general SU(N) case will be 



presented in |Tj]. Here we only quote the result for SU(2) 



, ^ 2 /4m 5 \ V6 5 /4m 5 \ 1/3 ^ , /2n , . 

r J {m) = l--l—j -gi^— J +0(m 5 / 2 ) (3.8) 

The numerical results for 77(771) are in the excellent agreement with this formula. 

Consider now the lower branch. i?i ow (0) and, due to Eq.( |3.5|) , also pi ow (0) go 
to zero in the limit m — > 0. Thereby, the solution looks very much similar to a 
combination of two real walls separated at a large distance (The phases a(z) and 
(3(z) are changed, of course, so that the whole solution is complex, but the change 
occurs in the central region where the absolute values of the fields R and p are 
very small.). One can say that it is a very loosely bound state of two real walls. 
The profile p(z) for the lower BPS branch at m = .5 found numerically is shown in 
Fig. ^[ Indeed, the complex wall is split apart in two widely separated real walls. 
p(0) ~ 6.7 • 10~ 4 in this case. We observe numerically that, for small masses, p(0) 
and R(0) go down oc exp{— C/m a } with a w 1/4, C p w 6.1, Cr w 5.3. 
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Figure 3: Upper BPS branch: p(z) for different masses. 




Figure 4: Upper BPS branch: R(z) for different masses. 



fO 




Figure 5: Profile p(z) of the lower BPS branch at m = 0.5 



The revealed phenomenon — the appearance of two solutions of the equation 
system (|3.2|) at m < m» and the absence of solutions at m > m* can be understood 
in the framework of the qualitative theory of differential equations (see Ref. || for a 
related discussion). The phase space of the system (|3.2|) is 4-dimensional. The space 
of all possible integral curves is 3-dimensional (there are 4 integration constants, but 
one of them corresponds to shifting the independent variable z along the trajectory). 

The system involves the integral of motion Q3.4|) so that the whole 3-dimensional 
space of the trajectories presents a set of 2-dimensional slices with a given value of 
Im W. 

We are interested in a special trajectory belonging to the slice with Im W = 



and passing through two chirally asymmetric minima (|1.9|) . These minima present 
singular points of the equation system ( |3.2| ) where the right hand sides of the equa- 
tions turn to zero. Thereby, the alias for our domain walls is separatrices. 

Linearizing the equations at the vicinity of singular points and studying the ma- 
trix of the linear system thus obtained (the Jacobi matrix), one can be convinced 
that these points have saddle nature: the Jacobi matrix has two positive A, \x and two 
negative —A, — /x eigenvalues. The space of all possible trajectories which have their 
origin at one of the singular points (say, the point with real 0, x) is 1 -dimensional 
- only the directions in the phase space along the eigenvectors of the Jacobi ma- 
trix with positive eigenvalues should be taken into account. There are two such 
eigenvectors which can be mixed with 1 real parameter. Likewise, the subspace of 
the trajectories which end up at the point <p — e -Mr//3 (3m/4) 1//6 , \ = ^(4/3m) 1 / 4 is 
also 1 -dimensional. They correspond to a mixture of two eigenvectors of the cor- 
responding Jacobi matrix with negative eigenvalues. The wall solution belongs to 
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a) m< m 



b) m = m 



c) m > m 



Figure 6: A space of trajectories with Im W = for different masses. The solid 
lines describe the subspace of separatrices originating at one of the minima while the 
dashed lines describe the subspace of separatrices ending up at another minimum. 

both subspaces. 

Generally, two 1-dimensional lines embedded in the 2-dimensional space (the 
space of all trajectories with Im W = 0) can cross at a finite number of points. 
What is this number and whether the lines cross at all depends on dynamics. We 
have found the presence of two intersections at m < m*, of one intersection at 
m = m* and the absence of intersections at m > m*. This picture illustrated 
schematically in Fig. || looks quite natural. Q 



8 It is instructive to compare the situation with that in the model analyzed recently in Ref. |l5| 
where a family of degenerate BPS wall solutions was observed. The model involved two chiral 
superfields with the superpotential 



BPS system involves here 4 singular points (4 vacua of the theory): two saddles, a stable and an 
unstable focuses. The BPS walls run from the instable focus to the stable one. Again, the families 
of the relevant separatrices are 1-dimcnsional. As in our case, the family of all trajectories is 
3-dimensional. Here, however, the wall solutions lie in the one-dimensional subspace of the real 
trajectories satisfying the two conditions Im <f> = and Im X = (all the vacua are real here). 
Hence, the space of the BPS wall solutions is 1 + 1 — (3 — 2) = 1-dimensional. 



2 

m 

W = — $ - 
A 
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Figure 7: Stationary points of the energy functional, m < m*. 

4 Wallsome Sphalerons. 

The BPS solutions realize the absolute minimum of the energy functional for all 
wall configurations. The other configurations have larger energy. In the region 
m < m*, we have found two different BPS solutions which should be separated by 
an energy barrier. Besides two nontrivial complex wall solutions, there is always a 
configuration describing two real walls separated at infinite distance. The energy 
of such a configuration is twice the energy of the real wall which coincides in our 
[S77(2)] case with the energy of the complex wall. There should be also an energy 
barrier separating this configuration from the lower BPS solution. This situation 
is illustrated in Fig. |j] where the energy functional minimized over all symmetric 
wall configurations with a fixed value of -R(O) (the value of \<f>\ in the middle) is 
drawn schematically as a function of R(0). We see that besides the minima (the 
BPS solutions), this plot has local maxima. In fact, they are the saddle points of 
the energy functional with only one unstable mode corresponding to changing R(0). 
Such saddle point configurations (the sphalerons) present non-trivial solutions of the 
equations of motion. They are not BPS solutions, of course. 

Thus, to study sphalerons, we should solve the full-scale second order equations 
of motion for the effective lagrangian with the potential ( |1.7| ) f\ and the standard 
kinetic term. Introducing again the polar variables x = pe ia , 4> = Re l/3 , the 

9 As usual, the sign of the potential should be reversed: we are looking for the minimum of the 
energy functional, not the action functional. 
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equations can be written in the form 



1 fi/? 5 

R"-R{3' 2 = 8R 3 \L{L + 3/2)+(3 2 + ] + — ^ - 4mi? 2 cos(/3_) 

3p z 

R(3" + 2R'(3' = l2R 3 (3 + + 4mR 2 sm(f3_) 

„ , 2 8R 4 2 16R G 

p - pa = L + m p — — 

p 9p 6 

8R\ 8mR 3 ,„ . , , 

pa" + 2p'a' = P+-— — sin(/3_) (4.1) 

p 3p 

where L = ln(i? 3 p 2 ), (3± = 3/3±2<x The system ( |4.1| ) involves one integral of motion 

T-U = R' 2 + R 2 (3' 2 + p' 2 + p 2 a' 2 - 

AR\L 2 + Pi ) + m 2 p 2 + ^ - ^ cos( ^_)j = const (4.2) 

9p z 3 

Solving the system (|4.1|) is a more complicated technical problem than that for the 
BPS system ( |3.2| ). The Cauchy problem involves here 8 initial conditions. As earlier, 
we can restrict our search with the symmetric configurations [see Eq.( j3.3|) ]. Then 4 
initial conditions are fixed: 

R'(0) = 0, p'(0) = 0, a(0) = tt/4, P(0) = -vr/6 



Four other initial values satisfy the relation (|4.2j ) with const = 0. Thus, we are left 
with 3 free parameters, say, p(0), R(0), and /3'(0), which should be fitted so that 
the solution approach the complex vacuum in Eq. fll.9p at z — > oo. Technically, we 
used the mismatch parameter in Eq. ( |3.6| ) as a criterium, the initial data were fitted 
to minimize it. 

The results are presented in Fig. We obtain first of all our old BPS solutions 
(solid lines in Fig.|8|). We find also two new solution branches drawn with the dashed 
lines in Fig. |8|. We see that, similarly to the BPS branches, two new dashed branches 
fuse together at some m = m** w 4.83. No solution for the system fl4.1|) exists at 
m > m**. For large masses, the minimal mismatch parameter A defined in Eq.( [3.6|) 
is not zero, and the larger is m — m**, the larger A is. Numerically, we find in this 
case 

Aeq.mot.M = 0.052(m - m„) a82 (4.3) 
in the full analogy with Eq.(| 



10 To be on the safe side, we would claim to** = 4.835 ± .020. The accuracy here is much lower 
than that for m» precisely because looking for the minimum of A in the space of 3 parameters 
is much more complicated numerical problem than fitting just 1 parameter which we did while 
searching for BPS solutions. That is also the reason why the dashed branches are not drawn in 
the region of small masses. Our accuracy is not sufficient there. 
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m 



Figure 8: The solutions of equations of motion at different masses. Solid lines 
correspond to the BPS solutions and dashed lines — to the sphalerons or non-BPS 
local minima. 

The physical interpretation of these new solutions is rather transparent. For 
m < m*, they present the sphalerons as shown in Fig. [?]. At m = m*, two BPS 
minima fuse together and the energy barrier separating them disappears. The upper 
sphaleron branch meets with the BPS branches at this point. When m is slightly 
above m*, the former BPS minimum is still a local minimum of the energy functional, 
but its energy is now slightly above the BPS bound (see Fig. ||a). The correspond- 
ing solution is described by the analytic continuation of the upper sphaleron branch. 
The lower dashed branch in the region m* < m < m** is still a sphaleron. At the 
second critical point m = m**, the picture is changed again (see Fig. ^b). The local 
maximum and the local minimum fuse together and the only one remaining station- 
ary point does not correspond to an extremum of the energy functional anymore. 
At larger masses, no non-trivial stationary points are left. 

Like in the simpler BPS case, the presence of the discrete number of the wall 
solutions for the system ( |4.1| ) depending on the value of the parameter m can be 
naturally interpreted in the language of the qualitative theory of differential equa- 
tions. The phase space is now 8-dimensional and the space of all trajectories is 
7-dimensional. There is one integral of motion (|4.2|) and our wall solutions belong 
to the 6-dimensional subspace of the trajectories with T — U = 0. Again, the 
minima Q1.9Q with zero derivatives = x = present the critical points of the 
equation system ([4.11) . The corresponding 8x8 Jacobi matrices involve 4 positive 
and 4 negative eigenvalues [the eigenvalues are actually the same as for the BPS 
system (|3.2|) , but each eigenvalue is now double degenerate]. Thereby, the family 
of separatrices which go down from one of the critical points is 3-dimensional. The 
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Figure 9: Stationary points of the energy functional: a) < m < m**, b) m = m**. 



same is true for the family of separatrices which end up at the other minimum. 
Generally, two 3-dimensional manifolds embedded in the 6-dimensional space of all 
trajectories with T — U — may cross at a finite number of points. We have seen 
that the crossing occurs at 4,2,1, or points depending on the value of mass. This 
looks reasonable though in this case we are not clever enough to visualize it. 

We hasten to comment that the height of the energy barrier as shown in Figs. 
[7|, [| is very much exaggerated. In fact, we could not even accurately determine 
the difference of the energy of the lower sphaleron branch and of all the dashed 
branches in the intermediate mass region m* < m < m** from the BPS bound epps- 
We can only say that the relative difference k = (e — £bps) / £bps is well below 10~ 4 
in this case. The appearance out of the blue of such a small number which is not 
present neither in supersymmetric QCD lagrangian nor in the effective TVY theory 
is rather surprising. Right now we do not see a profound theoretical reason for 
this "experimental fact" , but the numerical evidence is unequivocal. The fact that 
the barrier is very low conforms well with the fact that the second phase transition 
occurs pretty soon after the first one: the values m, ps 4.67 and m„ ps 4.83 are very 
close. 

The difference of the energy of the upper sphaleron from the BPS bound can 
be numerically evaluated in a certain mass region. The results are plotted in Fig. 
pL0| . Naturally, the height of the barrier goes to zero when m approaches m*. But, 
as seen from Fig. [10], even at small masses, the barrier height is still rather small. 
In the range of m* — m where we could calculate k reliably, it follows nicely the fit 
k « 3.90 • 10~ 5 exp{1.74(m„ - m)} (the dashed line in Fig. [TDD 

Finally, we present for illustration the profiles R(z) for all four branches at m = 3 
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Figure 10: Energy of the upper sphaleron as a function of mass. 



(see Fig. 0). 



5 Discussion. 



We have found that the properties of the system are drastically changed at m = 
4.67059 ... A and further at m ~ 4.83A. It makes sense to express the result in 
invariant terms and to trade A for an invariant physical quantity such as the gluino 
condensate S = | < Tr A 2 > | in a chirally asymmetric vacuum. From Eqs.(|1.9|), 



751) , we obtain 



167T 2 

^3 



Thus, the phase transitions occur at 



0.8£ 1/3 



(5.1) 



(5.2) 



Needless to say, they are not phase transitions of a habitual thermodynamic 
variety. In particular, the vacuum energy is zero both below and above the phase 
transition point — supersymmetry is never broken here. Hence E vac (m) = is not 
singular at m — m*, m**. Some similarities may be observed with the 2-dimensional 
Sine-Gordon model where the number of the states in the spectrum depends on the 
coupling constant (3 so that the states appear or disappear at some critical values of 
P ||16|| . May be a more close analogy can be drawn with the N = 2 supersymmetric 
Yang-Mills theory. The spectrum of that system depends on the Higgs expectation 



17 



z 



Figure 11: Profiles R(z) for BPS walls (solid lines) and wallsome sphalerons (dashed 
lines) for m = 3. 



value u = < Tr tp 2 >. A study of the exact solution of the model due to Seiberg and 



Witten [I7|| displays the existence of a "marginal stability curve" in the complex u- 
plane |fL8| . When crossing this curve, the spectrum pattern is qualitatively changed. 

The particular numerical value ( |5.2| ) was obtained by studying the effective TVY 
lagrangian with the superpotential ( |1.6| ) and the standard kinetic term \d(f)\ 2 + |<9x| 2 
which follows from the term | J d 4 #($$ + XX) in the superlagrangian. However, 
in contrast to the superpotential, the form of the kinetic term in the effective la- 
grangian cannot be fixed from the symmetry considerations alone. As as obvious 
generalization, one can write 



C 



kin 



a r 
4 J 



d 4 e<$>§ + ? / <fexx 



(5.3) 



Adding the term | / W($, X)d 2 + H.c. and excluding the auxiliary fields, the 
bosonic part of the lagrangian acquires the form 



C 



a\d^\ 2 + P\d X \' 



1 


dW 


2 1 


dW 


a 




P 


dx 



The BPS equations are modified to 

dd> . dW 



a 



dz 



dcj) 



az 



(5.4) 



(5.5) 



Suppose first that a = (3. Then the equations ( |5.5| ) can be reduced in the standard 
form (ll.lOQ by rescaling the independent variable z = az'. The walls become thicker 
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(or thinner) but, in all other respects, the form and the properties of the solutions 
are the same. If a ^ /3, the situation is not so trivial. The corresponding effective 
lagrangian seems to describe a different dynamics than that with a = (3 = 1. 



Nevertheless, we will shortly see that the lagrangian (^4) can still be reduced to the 
lagrangian with the standard kinetic term by rescaling both the independent variable 
z and the mass parameter m. 

Suppose, we have already set a = 1 by rescaling z in a proper way. It is not 
difficult to see that the BPS system (|5.5| ) is reduced to the standard system ( |1.10| ) 
with the superpotential 



2 ~ 
W = -$ 3 
3 



(5 - 6) 



where X' = (3 l ' 2 X. The expression ( |5.6| ) has the same functional form as Eq.( |1.6| ), 
only the scales for A and m have changed. Again, the BPS solutions disappear when 
m > m*(/3) = 4.67059. . ./5 6 / 5 . In physical units, m*(/3) « 0.8/3S 1 / 3 . The second 
critical point scales in the same way. 

Thus, we cannot claim that the phase transition in the theory of interest ( |1 . 1|) 
would occurs exactly at the mass ( |5.2j ). Our guess, however, is that the true values 
of m* and m** in the supersymmetric QCD are close to the estimate (|5.2j ). 

To summarize, our main physical conclusion is the absence of complex domain 
walls connecting different chirally asymmetric vacua in the supersymmetric QCD at 
large enough masses. In particular, they are absent in the limit m — > oo, for the 
pure supersymmetric Yang-Mills theory. 

In recent |Tj|, the complex domain walls in the pure N — 1 SYM theory were 
discussed in the context of D-brane dynamics. It was shown, in particular, that 
these walls, in spite of the fact that, on the fundamental level, they are composed of 
the fundamental gluon and gluino fields in the adjoint representation of the gauge 
group, behave like objects in the fundamental representation in color. In particular, 
a QCD string originating at a heavy fundamental source may end up at the domain 
wall (cf. the analogous situation in the massless adjoint QCD2 where fundamental 
colored sources are screened in spite of the fact that the lagrangian of the theory 



does not involve the fields which are fundamental in color f2(J). But that was in the 
assumption that such walls are there which, as we have seen, is not true. 

The walls interpolating between a chirally asymmetric and the chirally symmetric 
vacua are present at any value of mass and also in the limit m — > 00. For the time 
being, it is difficult for us to say something about their color properties: their very 



n Of course, the modification (5.3) is not the only possible one. One could invent more compli- 
cated kinetic terms involving higher (covariant) derivatives and some fractional powers of $ and 
X. We did not study these exotic modifications and cannot exclude a logical possibility that, after 
adding the terms with high derivatives, physics would dramatically change: e.g. the complex BPS 
walls would exist at any value of mass or would not exist at all. We do not believe, however, that 
this possibility is realized. 
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existence has been detected only in the framework of the effective TVY lagrangian 
involving only colorless composite fields. 

Finally, in the light of our new results, let us discuss the problem of torons. 
Torons are Euclidean gauge field configurations with fractional topological charge 
( |2.3|) , v being the integer multiple of 1/N, which may in principle contribute in the 
Euclidean path integral in the theory involving only adjoint color fields regularized 
in the infrared by placing it in a finite 4-dimensional toroidal box Whether 
such toron contributions are relevant or not in the supersymmetric gluodynamics is 



an old and controversial question (see e.g. |2^| ). It was discussed in details recently 
in [|J, and we will not present here anew all the arguments pro and contra. 

We would like only to reiterate here that the presence of torons would substan- 
tially modify the physical picture for the vacuum structure of the supersymmetric 
gluodynamics (vis. N + 1 degenerate vacua and N different real domain walls in- 
terpolating between them) advocated in M and in the present paper. In the first 
place, the presence of torons would modify the gluing prescription for the effective 
potential (|2.2| ). The gluing would occur at arg(0 3 ) = 7r rather than at arg(0 3 ) = n/2 
and arg(0 3 ) = 37r/2 as is the case when only the instanton-like gauge field configu- 
rations with integer topological charge are taken into account. Then the potential 
( |2.2| ) would have only one non-trivial minimum at arg(0 3 ) = 0, not two such minima 
as before. The theory would involve only two vacua, the chirally asymmetric and 
the chirally symmetric one, and only one interpolating domain wall. 

We have seen, however, that, at least if the supersymmetric gluodynamics is 
understood as the supersymmetric QCD in the limit m — ► oo, this is not true. 
When the matter chiral multiplets are added, torons just do not exist however large 
the mass is (the Dirac equation for fundamental color fields does not make sense on 
a toron background), and the picture with 3 [N + 1 for the SU(N) case] vacuum 
states is correct. 

A nontrivial feature of the theory under consideration is that the real domain 
walls do not decouple in the limit m — ► oo. That is in contrast to the simple 
two-dimensional model discussed in Ref. M where the energy of a two-dimensional 
analogue of the 4-dimensional supersymmetric domain walls becomes infinite in the 
limit m — > oo so that, in this limit, the vacua decouple completely from each other, 
and the physical picture corresponds not to the spontaneous breaking of a discrete 
symmetry, but to the appearance of a new supers election rule in the theory which 
just corresponds to the analogue of toron contributions in the path integral being 
switched on. 

Note further that, for the supersymmetric QCD with large mass, the problem 
of how the different branches of logarithm in the potential are glued together is 
completely irrelevant as far as the wall solutions are concerned. For the real walls, 
we have always arg(0 3 % 2 ) = irrespectively of whether arg(0 3 ) = (for the wall 
connecting the vacuum with the positive value of < Tr A 2 > with the chirally 
symmetric one) or arg(0 3 ) = it ( for the wall connecting the vacua with < Tr A 2 > 
< and < Tr A 2 > = 0). The phase of the heavy field x always compensates the 
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phase of 0. In a sense, heavy fields do not decouple completely in this case. 

The question of whether it still makes sense to formulate the pure super symmetric 
gluodynamics not as the large mass limit of supersymmetric QCD but in some other 
way which takes into consideration the toron contributions so that the theory would 
involve only one domain wall connecting two different vacua is still open. To clarify 
it, one should first try to understand better the dynamics of the chirally symmetric 
phase. 

In fact, the existence of the chirally symmetric vacuum is probably the most 
interesting and unexpected result of the effective lagrangian approach. What is 
especially surprising is that it has the same (zero) energy as the chirally asymmetric 
vacua. Such degeneracy is not required by any symmetry considerations. Note that 
the degeneracy is there only for the theory defined at infinite volume. At small 
volumes when the effective gauge coupling is small, the theory can be analyzed 
directly in terms of the fundamental gluon and gluino fields, there are only N vacuum 
states with nonzero value of the condensate, and no trace of the chirally symmetric 
phase is seen [p3| . By analyticity, that should be true also for large (but finite) 
volumes. Indeed, playing around with the finite volume formulas of Ref. displays 
that the degeneracy is lifted also for large volumes, the energy gap going to zero 
exponentially fast in the limit L — > oo. 

There is a striking analogy with a two-dimensional model [QCD2 with the 
fermions in the adjoint color representation for higher unitary groups SU(N > 3]) 
studied in Ref.||13||. In this model, the bosonized lagrangian displays the existence 
of several vacuum states separated presumably by the domain walls (the latter are 
just solitons in two dimensions). Like in the supersymmetric QCD4, different vacua 
are distinguished by the different values of the fermion condensate < Tr {XrXl} >■ 
For higher even N, the existence of several vacua does not follow from symmetry 



considerations (for odd N, it does). That was formulated as a paradox in Ref. fl3 
but, bearing in mind fresh insights from supersymmetric QCD^, it looks now less 
surprising. 
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